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摘要 : 为 探究 Astrometry.net 和 SCAMP 两 款 软 件 在 天 测 位 置 定 标 时 的 应 用 差异 以 及 以 不 
同方 式 运行 软件 所 得 定 标 效 果 的 差异 ， 分 别 介 绍 了 两 款 软件 位 置 定 标的 原理 ， 以 ZTEF(Zwicky 
Transient Facility) 的 100 张 巡 天 图 像 为 测试 图 像 ， 以 Gaia DR2 为 参考 星 表 ， 分 别 以 4 种 不 同 
的 定 标 流程 对 测试 图 像 进行 了 位 置 定 标 。 从 计算 总 耗 时 、 星 表 匹 配 情况 和 坐标 均 方 根 ( Root 
Mean Square，RMS ) 偏差 3 个 方面 ， 比 较 了 不 同 流程 所 得 定 标的 结果 ， 同 时 也 与 ZTF 给 出 的 
定 标 结果 进行 了 对 比 。 结 果 表 明 Astrometry.net 能 够 快速 给 出 粗略 的 线性 定 标 结果 ，SCAMP 
可 以 快速 在 前 者 的 基础 上 直接 进行 扭曲 修正 ， 两 者 相 结 合 的 定 标 流程 平均 耗 时 仅 为 1s， 所 
得 定 标 结果 的 平均 赤 经 赤 纬 均 方 根 偏 差 均 小 于 70 毫 角 秒 (小 于 0.1 个 像素 ) ， 优 于 仅 使 用 
Astrometry.net 所 得 定 标 结果 。 由 此 可 见 ， 搭 配 使 用 Astrometry.net fe SCAMP 可 以 实现 更 快速 
且 更 准确 的 天 测 位 置 定 标 ， 可 推广 应 用 到 时 域 天 文 巡 天 项 目的 数据 处 理 上 。 
关键 词 : 天 体 测 量 软件 研究 ; CCD 图 像 处 理 ; Astrometry.net; SCAMP 
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光学 波段 天 文 图 像 通 和 常 由 望远镜 配备 的 电子 厢 合 絮 件 (Charge Coupled Device, CCD) 芯片 相机 采 
集 ， 主 要 采用 FITS( Flexible Image Transport System) 格式 进行 存储 。 为 了 测量 图 像 上 特定 星 像 的 亮度 
或 视 星 等 ， 我 们 一 般 采 用 较 差 测 光 或 图 像 相 减 的 方式 一 。 对 于 前 者 ， 先 在 观测 图 像 ( 测 光 图 像 ) 上 选 
取 一 些 参考 星 ， 通 过 比较 参考 星 的 仪 硕 星 等 和 这 些 参 考 星 在 这 天星 表 中 的 视 星 等 ， 得 到 图 像 零点 ， 有 再 
对 上 述 特定 源 的 仪 硕 星 等 进行 零点 修正 ， 从 而 得 到 该 源 基 于 巡天 星 表 的 视 星 等 。 对 于 选取 的 参考 星 ， 
我 们 需要 通过 其 在 CCD 网 像 中 的 像素 坐标 计算 出 对 应 的 天 球 坐 标 ， 再 与 巡天 星 表 中 的 对 应 目标 进行 
匹配 。 观 测 图 像 的 位 置 定 标 越 准 确 ， 图 像 中 星 像 与 参考 星 表 中 恒星 的 坐标 匹配 结果 越 可 牧 ， 否 则 可 能 
出 现 误 匹配 (图 像 某 星 像 与 匹配 上 的 星 表 中 的 恒星 不 是 同一 颗 星 ) 的 情况 ， 从 而 影响 零点 的 准确 度 。 
同时 ， 获 得 并 报告 特定 源 的 天 球 坐 标 信息 也 是 测 光 数据 处 理 的 主要 目标 之 一 。 显 然 ， 观 测 图 像 上 各 是 
像 的 位 置 定 标 是 整个 测 光 数据 处 理 过 程 中 非常 重要 的 一 环 。 

从 图 像 像 素 坐 标 到 天 球 坐 标的 转换 关系 的 系数 以 及 相关 信息 记录 为 FITS 文件 头 中 的 一 系列 特定 
关键 字 的 数值 和 备注 信息 ， 这 些 关 键 字 称 为 图 像 的 世界 坐标 系 | ( World Coordinate System, WCS) 信 
上 县。 根据 文 L4] ， 从 图 像 像 素 坐 标 到 天 球 坐 标的 转换 可 以 分 为 3 步 : (1) 从 像素 坐标 (Pixel Coordinates) 
转换 到 以 度 为 单位 的 投影 平面 坐标 (Projection Plane Coordinates); (2) 从 投影 平面 坐标 转换 到 本 地 球 
面 坐标 (Native Spherical Coordinates); (3) 从 本 地 球面 坐标 转换 到 最 终 的 天 球 坐 标 ( Celestial Spherical 
Coordinates) 。 第 2 步 操 作 视 投影 类 型 而 定 ， 有 固定 的 转换 公式 ,在 图 像 中 常 使 用 心 册 正切 投影 
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(TAN ) 。 第 3 步 为 不 同 球面 坐标 系 下 的 坐标 转换 ， 转 换 公式 也 是 固定 的 。 因 此 ， 天 测 位 置 定 标 实 
际 要 计算 的 就 是 第 1 步 中 的 转换 系数 。 

理想 情况 下 ， 从 像素 坐标 到 投影 平面 坐标 的 转换 可 以 由 线性 变换 近似 描述 ， 包 括 坐 标的 平移 、 旋 
转 、 伸 缩 、 切 变 等 ,但 由 于 大 气 和 仪器 等 因素 (如 大 气 折射 ”、 大 气 淇 流 、 望 远 镜 CCD 缺陷 、 光 路 缺 
陷 和 重力 形变 等 ) ， 望 远 镜 拍摄 到 的 图 像 相 比 理想 的 投影 平面 存在 一 定 程度 的 扭曲 。 为 了 保证 图 像 各 
个 区 域 坐标 具有 足够 的 准确 性 ， 一 般 的 方法 是 引入 高 阶 多 项 式 作 为 修正 项 拟 合 图 像 扭曲 ， 以 修正 坐标 
偶 差 。 帮 望远镜 视 场 较 小 且 仪 硕 状 态 恨 好 ， 则 对 所 担 摄 图 像 定 标 时 仅 计 算 线 性 转换 系数 (后 文 简称 线 
性 定 标 ) 就 可 以 满足 需求 ， 但 对 于 巡天 图 像 这 类 视 场 较 大 的 图 像 ， 边 绿 部 分 扭曲 比较 明显 ， 此 时 就 需 
要 额外 计算 高 阶 修正 系数 (后 文 简称 扭曲 修正 ) 以 得 到 理想 的 定 标 结 采 。 

当前 ，Astrometry.net “| 和 SCAMP 是 使 用 较 多 的 两 款 天 测 位 置 定 标 软件 ， 前 者 主要 特点 为 在 
没有 任何 关于 图 像 位 置信 息 时 可 以 计算 一 套 完整 的 世界 坐标 系 信息 ， 而 后 者 主要 功能 为 对 一 套 初 始 
的 世界 坐标 系 信息 进行 修正 。 本 文 首先 介绍 两 蒜 软 件 的 定 标 原理 ， 接 着 以 ZTF( 美 国 的 一 项 正在 进 
行 中 的 时 域 巡天 项 目 ) 巡 天 ”的 科学 图 像 为 测试 图 像 ， 计 算 并 比 对 了 在 输入 源 表 相 同 的 情况 下 ， 以 4 
种 不 同方 式 运 行 两 款 软 件 得 到 的 天 测定 标 质量 。 文 中 使 用 的 Astrometry.net 版 本 为 0.89，SCAMP 版 本 
为 2.10. 


1 天 测 位 置 定 标 原理 介绍 


如 前 所 述 ， 从 像 系 坐 标 到 天 球 坐 标 转换 的 关键 在 于 像素 平面 坐标 到 投影 平面 坐标 的 转换 。 硅 不 考 
虑 扭曲 修正 ， 则 两 者 之 间 的 转换 关系 由 FITS 头 文件 中 的 CD, ,给 定 ,i 和 j 取 值 1 或 2。 线 性 转换 公式 为 
x CD,, CD,, u 
CI 
CD,, CD,, v 
RP, (u, 0) 为 以 参考 点 为 原点 的 像素 坐标 ; (x, y) 为 以 同一 参考 点 为 原点 的 投影 平面 坐标 ， 参 考点 
由 关键 字 CRPIX1 和 CRPIX2 给 定 。 大 考虑 扭曲 修正 ， 依 据 所 使 用 的 扭曲 修正 类 型 ， 先 修正 像素 坐标 
(u, v) 再 进行 线性 转换 ， 或 者 先进 行 线性 转换 再 修正 投影 平面 坐标 (*, y) 。 

天 测定 标 软件 一 般 通 过 图 像 星 表 和 参考 星 表 之 间 的 比较 分 析 来 计算 坐标 转换 中 的 各 项 系数 。 对 于 
Astrometry.net， 使 用 者 既 可 以 输入 图 像 由 软件 完成 找 星 生成 图 像 星 表 ， 也 可 以 直接 输入 准备 好 的 图 像 
星 表 由 软件 直接 进行 定 标 。 而 对 于 SCAMP ， 使 用 者 需要 提前 使 用 软件 SExtractor ”生成 图 像 星 表 ， 再 
输入 SCAMP 进行 定 标 修正 。 

1.1 Astrometry.net 

Astrometry.net 是 用 C 语言 编写 的 开源 免费 软件 ， 可 以 在 Linux 系统 和 MacOS 系统 上 安装 ( 见 
https ://astrometrynet.readthedocs.io/en/latest/) ， 安 装 后 可 以 通过 solve-field 命令 运行 。 该 软件 通过 匹 
配 图 像 星 表 中 的 星 标 ( asterism ) 和 参考 星 表 中 的 星 标 计算 图 像 的 世界 坐标 系 信息 。 一 个 星 标 由 4 颗 星 
组 成 ， 根 据 这 4 颗 星 的 坐标 可 以 计算 出 一 个 代表 它们 之 间 相 对 位 置 关 系 的 喻 希 人 码 ; 如 果 这 4 颗 星 既 在 
图 像 星 表 中 也 在 参考 星 表 中 ， 则 由 它们 的 像素 坐标 计算 的 哈 希 码 和 由 他 们 在 参考 星 表 中 的 天 球 坐 标 计 
算 的 哈 硕 人 码 是 接近 的 (由 于 存在 位 置 测 量 误差 和 扭曲 ， 它 们 的 哈 希 人 码 不 完全 相同 )。 

Astrometry.net 根据 不 同 的 光学 了 巡天 星 表 构 建 不 同 尺度 的 遍布 全 天 的 星 标 ， 这 些 星 标的 哈 希 码 以 及 
组 成 相应 星 标的 星 的 天 球 坐 标 存储 为 一 系列 index 文件 (http://data.astrometry.net/ ) 以 供 查 找 比 对 。 
Astrometry.net 根据 亮度 排序 ) 依次 对 图 像 星 表 中 的 星 像 构建 星 示 并 计算 哈 希 人 码 在 index 文件 中 寻找 
哈 希 人 码 相 近 的 星 标 作为 匹配 星 标 。 Astrometry.net 根据 找 到 | 的 匹配 星 标 计 算 世 界 坐标 系 AAR ， 据 此 计算 
图 像 上 其 他 星 的 天 球 坐 标 ， 并 与 index 文件 中 其 他 星 的 天 球 坐 标 进 行 比 较 ， 通 过 贝 叶 斯 决断 方法 判断 
该 匹配 是 否 为 真 。 夺 匹配 为 真 ， 则 完成 计算 ; 寿 为 假 ， 则 继续 验证 其 他 的 匹配 星 标 。 

作为 可 选项 ，Astrometry.net 可 以 根据 匹配 上 的 星 计算 图 像 扭曲 修正 系数 ， 采 用 的 扭曲 修正 多 项 式 
H SIP (Simple Imaging Polynomial) "™ ， 修 正 对 象 为 像素 坐标 (ww， v), ATA 
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C) CD, , CD, , c + flu, D) 
= x f (2) 
CD,, CD,, v +g(u, v) 
HP, flu, v) Al g(u, wv) 分 别 为 两 方 品 坐标 的 修正 多 项 式 函 数 ， 
flu, wD... Apat P + q SA opp, (3) 
g(u, v) =); 本 Pip gl» P +q SB orr, (4) 


A o ll B win 为 两 个 方向 的 修正 阶 数 ， 一 般 是 相等 的 ， 可 以 通过 命令 行 参 数 设置 A, MB, We 
项 式 中 的 系数 。 
1.2 SCAMP 

SCAMP 是 用 C 语言 编写 的 开源 免费 软件 ， 可 以 在 Linux 系统 和 MacOS RA EZZ (IL https: // 
github.com/astromatic/scamp/) ， 安 装 后 可 以 通过 scamp 命令 运行 。 该 软件 需要 使 用 者 提供 一 个 初始 的 
地界 坐标 系 信 息 ， 通 过 独立 的 两 步 完 成 对 世界 坐标 系 信 息 的 修正 。(1) 模 式 匹 配 (Pattern matching) , 
目的 是 对 初始 世界 坐标 系 中 的 线性 定 标 系 数 进行 修正 ; (2) 计 算 天 测 解 (Astrometric solution) ， 目 的 是 
在 给 定 线性 定 标的 基础 上 计算 扭曲 修正 。 使 用 者 可 以 根据 需要 ， 只 运行 其 中 一 步 ， 或 者 两 步 都 运行 。 

在 模式 匹配 中 ，SCAMP 首先 依据 初始 世界 坐标 系 信 息 ， 把 参考 星 表 中 的 天 球 坐 标 投影 转换 为 像 
素 坐 标 。 通 过 互相 关 分 析 比 较 图 像 星 表 像 素 坐标 对 和 参考 星 表 像素 坐标 对 的 间距 和 方位 角 的 二 维 分 布 
( 见 文 [7] 第 3 市 )，SCAMP 可 以 得 到 初始 世界 坐标 系 中 的 像素 比例 尺 和 旋转 角度 的 修正 量 ， 并 据 此 
更 新 世界 坐标 系 和 参考 星 表 像素 坐标 。 通 过 同样 的 方法 比较 图 像 星 表 像 素 坐 标 和 参考 星 表 像素 坐标 在 
像素 平面 上 的 分 布 ，SCAMP 可 以 得 到 两 套 坐 标 之 间 的 平移 偏差 ， 从 而 对 初始 世界 坐标 系 中 的 指 问 信 
县 进行 修正 。 

在 计算 天 测 解 中 ，SCAMP 对 图 像 星 表 和 参考 星 表 的 源 进行 匹配 ， 通 过 最 小 化 加 权 坐 标 差 值 平方 
和 ( 见 文 [7] 第 4 节 ) 的 方法 计算 扭曲 修正 系数 。SCAMP 采用 的 扭曲 修正 多 项 式 为 TPV, BEN% 
最 高 可 设 为 7 阶 ， 修 正 对 象 为 投影 平面 坐标 (x*, y) ， 修 正 公 式 见 文 [5] 中 的 (1) 式 和 (2) 式 。 


2 READ AAT LES 


2. 1 测试 图 像 选择 

本 文 随机 选择 了 100 个 均匀 分 布 于 北 天 的 坐标 ， 在 每 个 坐标 位 置 处 选择 一 张 具 有 较 好 视 宁 度 和 较 
深 极限 星 等 的 rz 波段 ZTF 巡天 科学 图 像 作 为 测试 图 像 (下 载 方式 见 https ://irsa. ipac. caltech. edu/docs/ 
program_interface/ztf_api.html ) 。 每 张 图 像 的 视 场 约 为 52’x52’, RRR TA 3072x3080, 像素 比例 尺 
为 1.012 角 秒 / 像 素 。 图 1 展示 了 测试 图 像 的 半 高 全 宽 、 极 限 星 等 (与 参考 星 表 较 差 后 对 应 波段 的 视 星 
等 ) 以 及 用 于 天 测定 标的 源 ( 见 2.2 节 ) 数 量 分 布 直方 图 。 半 高 全 宽 和 极限 星 等 均 从 图 像 的 FITS 头 文 
件 中 获取 。 

每 张 ZTF 科学 图 像 的 FITS 头 文件 中 市 有 准确 的 世界 坐标 系 信息 。 为 了 量化 测试 图 像 的 扭曲 程度 ， 
我 们 利用 图 像 的 世界 坐标 系 信息 ， 计 算 每 个 像素 坐标 考虑 扭曲 修正 和 不 考虑 扭曲 修正 分 别 对 应 的 天 球 
坐标 之 间 的 偏差 ， 取 所 有 像素 位 置 的 偏差 距离 的 最 大 值 作 为 图 像 扭 曲 程度 的 量化 指标 。 图 2 展示 了 随 
机 选取 的 一 张 图 像 的 扭曲 情况 。 在 此 指标 下 ， 所 有 测试 图 像 中 扭曲 程度 最 大 为 1. 46"( 约 1. 44 像素 ) ， 
平均 为 0.75”( 约 0.74 像素 )。 图 3 展示 了 图 像 扭 曲 程度 的 分 布 直方 图 。 

2.2 定 标 流程 和 参数 设置 

首先 删除 图 像 已 有 的 世界 坐标 系 信息 ， 使 用 SExtractor 对 图 像 进 行 找 源 ， 并 基于 该 软件 给 出 的 每 
个 源 的 测量 结果 筛选 适合 用 于 定 标的 源 。 找 源 时 的 重要 参数 设置 见 表 1。 筛选 时 ， 对 源 的 半 高 全 宽 进 
ÍT 3 FRE ZS Hi eV RAB a ST OR UPR BE FLAGS 为 0 的 源 以 去 挥 测量 结果 可 能 不 准确 的 
源 ， 同 时 去 掉 椭 率 大 于 0.5、 信 噪 比 低 于 10 的 源 ， 最 后 得 到 用 于 天 测定 标的 图 像 星 表 。 
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direction of distortion 


#1 SExtractor 重要 参数 设置 


Table 1 The configuration of major parameters for SExtractor 


Parameter Value Explanation 
DETECT_MINAREA 5 The minimum area above the detection threshold is required 
DETECT_THRESH 5.0 Detection target threshold 
ANALYSIS_THRESH 5.0 Analyze the target threshold 
FILTER N Use filter or not 
PHOT_AUTOPARAMS Zed 3:3 MAG AUTO parameter, Kron factor and minimum radius 


依次 通过 4 种 不 同 流程 进行 位 置 定 标 ， 分 别 为 : 

(1) 运 行 solve-field 命令 ， 一 次 性 进行 线性 定 标 和 扭曲 修正 ; 

(2) 首先 运行 solve-field 进行 线性 定 标 ( 使 用 -T 参数 ) ， 再 次 运行 solve-field， 在 第 1 次 运行 的 基础 
上 (使 用 -v 参数 ) 进行 包括 扭曲 修正 在 内 的 计算 ; 
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(3) 首先 运行 solve-field 进 WA 接着 运行 SCAMP 直接 进行 扭 遇 修 正 ; 

(4) 首 先 运 行 solve-field 进行 线性 定 标 ， 接 着 运行 SCAMP 进行 线性 定 标的 修正 ， 最 后 运行 
SCAMP 进行 扭曲 修正 。 

计算 的 扭曲 修正 均 为 4 阶 ， 与 图 像 自 带 的 修正 阶 数 一 致 。 这 4 种 流程 按 次 序 分 别 编号 为 流程 1 到 
流程 4。 不 同 流程 中 运行 solve-field 时 固定 的 命令 行 参数 选取 和 设置 见 表 2。 不 固定 的 参数 只 有 搜索 中 
心 坐 标的 设置 、 是 否 计算 扭曲 修正 及 其 阶 数 的 设置 和 是 否 验 证 已 有 世界 坐标 系 的 设置 。 不 同 流程 中 运 
ÍT SCAMP 时 国定 的 重要 参数 设置 见 表 3， 不 固定 的 只 有 是 否 进 行 模式 匹配 和 扭曲 修正 。Astrometry. 
net 的 参考 index 文件 选 为 基于 Tycho-2 星 表 和 Gaia DR2 星 表 的 5200 系列 ，SCAMP 使 用 的 参考 星 表 为 
Gaia DR2 星 表 。 所 有 计算 过 程 在 内 存 为 16 GB, 、 世 片 为 Apple M1 Pro 的 14 英寸 MacBook Pro 笔记 本 电 
脑 上 运行 。 


表 2 solve-field 命令 行 参数 设置 


Table 2 The configuration of solve-field command line options 


Command line parameter Value Explanation 
——scale-low 0.7 Lower limit of image field of view 
——scale-high 0.9 Upper limit of image field of view 
——scale-units degwidth The unit of the upper and lower limits of the field of view 
——radius 0.8 Search radius 
——crpix-center 一 一 Specify the center of the image as the reference point 
-一 width 3 072 Image width (pixel) 
——height 3 080 Image height (pixel) 
—-x-column XWIN_IMAGE Enter the column name for the horizontal coordinate in the image star list 
--y-column YWIN_IMAGE Enter the column name for the vertical coordinate in the image star list 
——sort-column FLUX_AUTO The name of the column to sort the image star list 


:该 命令 行 参 数 不 用 取 值 ， 选 定 即 可 。 


#3 SCAMP 重要 参数 设置 
Table 3 The configuration of major parameters of SCAMP 


Parameter Value Explanation 
ASTREFMAG_LIMITS 12.0, 18.0 Refer to the magnitude range used for scaling in the star list 
PIXSCALE_MAXERR 1.2 Maximum error multiple of pixel scale 
POSANGLE_MAXERR 5.0 Maximum azimuth deviation (degree ) 
POSITION_MAXERR 1.0 Maximum position deviation (Angle minute) 
PROJECTION_TYPE SAME Projection type 

DISTORT_KEYS XWIN_IMAGE, The column name in the image star list on which the distortion 
YWIN_IMAGE correction depends 
DISTORT_DEGREES 4 Distortion correction order 
NTHREADS 1 Number of running threads 


2.3 计算 定 标 质 

考虑 到 不 ras 环境 对 天 测定 标的 准确 性 要 求 不 同 ， 本 文 使 用 定性 和 定量 两 个 指标 衡量 天 测定 标 
的 质量 。 定 性 指标 选 为 定 标 后 图 像 星 表 源 坐标 是 否 在 全 图 范围 内 与 参考 星 表 源 坐 标 相 匹配 ， 匹 配 阔 值 
选 为 0.5 倍 图 像 半 高 全 宽 ， 用 于 匹配 的 参考 星 表 选 为 Gaia DR2 星 表 。 定 量 指标 则 选 为 定 标 后 的 图 像 
源 坐 标 与 参考 星 表 中 对 应 源 坐标 的 均 方 根 偏差 。 具 体 计 算 方法 如 下 。 
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(1) 使 用 SExtractor 对 定 标 后 图 像 分 别 以 高 低 国 值 进行 两 次 找 源 ， 高 探测 阀 值 和 分 析 姜 值 选 5. 0， 
低 探 测 装 值 和 分 析 闪 值 选 1.0。 以 与 2. 2 节 中 相同 的 筛选 条 件 对 高 浆 值 星 表 进行 澳 选 ， 再 与 参考 星 表 
进行 坐标 匹配 ， 接 着 将 参考 星 表 与 低 靖 值 图 像 星 表 匹 配 。 如 果 与 高 国 值 星 表 源 配对 的 参考 星 表 源 在 低 
国 值 星 表 中 配对 的 源 不 是 当前 高 国 值 星 表 源 ， 则 判定 为 误 匹 配 。 在 高 国 值 星 表 中 去 掉 所 有 误 匹 配 的 
源 ， 计 算 余 下 源 坐 标 与 匹配 上 的 参考 星 表 源 坐标 之 间 的 距离 ， 作 为 图 像 星 表 源 坐标 的 偏差 。 需 要 指出 
的 是 ，SExtractor 不 能 识别 由 SIP 多 项 式 描述 的 扭曲 修正 关键 字 ， 因 此 ， 对 于 Astrometry.net 给 出 的 定 
标 结 果 ， 我 们 先 使 用 python 包 sip_tpv( https ://github.com/stargaser/sip_tpv ) 将 SIP 关键 字 转 换 为 TPV 
KEF, PEH SExtractor 进行 找 源 。 

(2) 对 图 像 划 分 为 10 行 10 列 的 网 格 阵列 ， 将 每 个 网 格 内 源 的 坐标 偏差 除 以 0. 5 ARREARS, 
再 剔除 5 倍 标准 差 ( 标 准 差 选 为 84. 13 百 分 位 与 15. 86 百 分 位 差 值 的 一 半 ) 之 外 的 异常 值 后 ， 取 平均 值 
作为 当前 网 格 内 源 坐 标准 确 性 的 衡量 ， 最 后 得 到 的 数组 即 可 反映 图 像 各 区 域 坐标 的 偏差 情况 。 图 4 展 
示 了 一 张 图 像 只 进行 线性 定 标 和 扭曲 修正 后 的 坐标 偏差 图 。 如 果 数 组 内 所 有 值 小 于 1， 则 认为 图 像 星 
表 在 全 图 范围 内 与 参考 星 表 相 匹配 ， 定 性 为 定 标 质 量 好 ; 如 果 数 组 内 大 于 1 的 值 的 个 数 大 于 0 HAE 
过 10， 则 定性 为 定 标 质 量 中 等 ;否则 定性 为 定 标 质 量 差 。 
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图 4 图 像 各 区 域 坐标 偏差 图 。(a) 仅 线性 定 标 ; (b) 扭曲 修正 后 定 标 


Fig.4 Coordinate deviation map. (a) Linear calibration; (b) distortion corrected 
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(3) 计 算 剔除 异常 值 后 所 有 源 的 赤 经 赤 纬 坐标 偏差 ( 定 标 后 的 坐标 减 去 参考 星 表 的 坐标 ) 以 及 相应 


的 均 方 根 ， 作 为 该 图 像 天 测定 标 质量 的 定量 指标 。 图 5 展示 了 所 有 图 像 上 的 源 只 进行 线性 定 标 和 进行 
扭曲 修正 后 的 坐标 偏差 分 布 图 。 
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AS 所 有 图 像 源 赤 经 赤 纬 偏差 分 布 。(a) 仅 线性 定 标 ; (b) 捏 曲 修正 后 定 标 
Fig.5 2D-distributions of deviation along R.A. and Dec. (a) Linear calibration; (b) distortion corrected 


除了 以 上 提 到 的 指标 ， 赤 经 赤 纬 平均 偏差 及 偏差 的 标准 差 与 星 等 的 关系 图 也 能 直观 地 展示 定 标 质 


量 的 好 坏 。 图 6 展示 了 所 有 图 像 上 的 源 在 扭曲 修正 后 赤 经 赤 纬 仿 差 随 星 等 的 分 布 图 ， 包 括 各 星 等 处 的 
偏差 中 值 及 1 EAE TE, FA 6 可 见 ， 两 方 问 的 坐标 偶 差 在 各 星 等 处 较 均 义 地 分 布 在 零点 左右 。 
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图 6 所 有 图 像 源 偏 差 随 星 等 分 布 图 。 以 点 连接 的 实 线 代 表 不 同 星 等 处 的 坐标 偏差 中 值 ， 
阴影 部 分 代表 坐标 偏差 的 1 REA (a) A; (b) 杰 续 
Fig.6 Distribution map of all image source deviations with magnitude. The median deviation at different magnitude 


are plotted as solid lines with lo boundary marked as dashed lines. (a) R.A.; (b) Dec 
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3 定 标 结果 比 对 


本 市 从 总 计算 耗 时 、 定 性 的 星 表 匹 配 情况 和 定量 的 坐标 均 方 根 仿 差 3 个 方面 对 不 同 定 标 方式 得 到 
的 结果 进行 比较 。 比 较 对 象 除了 2. 2 市 中 提 及 的 4 种 定 标 流程 的 结果 ， 也 包括 Astrometry.net 仅 线 性 
定 标的 结果 和 SCAMP 仪 线 性 定 标 修正 的 结果 ， 以 及 ZTF 图 像 已 有 的 定 标 结果 。 表 4 记录 了 这 些 定 标 
结 末 的 平均 计算 耗 时 ， 定 标 质量 分 别 为 好 、 中 等 、 差 的 图 像 占 比 ， 以 及 100 张 图 像 平均 赤 经 赤 纬 方 癌 
的 均 方 根 偏差 。 


表 4 不 同 定 标 方式 的 耗 时 、 定 性 质量 统计 和 平均 赤 经 赤 纬 均 方 根 偏差 
Table 4 Time cost, quality statistic and average RMS deviation along R.A. and Dec axis of 


calibration results from different processes 


Average Fraction of Fraction of Fraction of Average RMS Average RMS 


Calibration process time cost good quality medium quality bad quality along R.A. along Dec 
/s /(%) /(%) /(%) /mas /mas 
T by 0.72 29 40 31 329. 18 314. 60 
cee re aes 4. 59 68 32 0 184. 17 169. 82 
Process 1 1. 64 83 8 9 177. 81 173. 26 
Process 2 5. 95 96 4 0 73. 39 69. 99 
Process 3 1. 00 100 0 0 68. 68 64. 36 
Process 4 5. 59 100 0 0 67. 80 63. 69 
ZTF calibration — 100 0 0 68. 66 64. 58 


3.1 线性 定 标 修正 前 后 结果 对 比 

定 标 流程 4 中 包含 SCAMP 对 Astrometry.net 给 出 的 线性 定 标 结果 进行 修正 的 步骤 ， 因 此 可 以 对 修 
正 前 后 的 线性 定 标 质量 进行 比较 。 图 7 展示 了 修正 前 后 定 标 质 量 分 别 为 好 、 中 等 和 差 的 图 像 数 量 分 
布 。 图 7 表明， 经 过 SCAMP 修正 后 ， 定 标 质量 好 的 图 像 占 比 由 29% 升 高 至 68% ， 而 质量 差 的 图 像 占 
比 由 31% 降 至 O%, Al 8 展示 了 修正 前 后 杰 经 杰 纬 方向 均 方 根 仿 差 的 分 布 。 由 图 8 可 以 看 出 ， 经 过 修 
正 后 ， 两 方向 均 方 根 俩 差分 布 更 集中 且 均 值 更 小 。 图 9 展示 了 运行 solve-field 进行 线性 定 标 和 单线 程 
运行 SCAMP 进行 线性 定 标 修 正 的 耗 时 分 布 。 由 网 9 可 以 看 出 ， Astrometry.net 一 般 在 1s 内 即 可 完成 
图 像 的 线性 定 标 ， 而 SCAMP 则 需要 4~5 s 才能 完成 线性 定 标的 修正 。 
3.2 4 种 流程 定 标 结果 对 比 

我 们 对 4 种 不 同 定 标 流程 所 得 结果 的 质量 进行 对 比 。 图 10 展示 了 4 种 流程 定 标 质 量 分 别 为 好 、 
中 等 和 差 的 图 像 数量 分 布 。 由 图 10 可 以 看 出 ,流程 1 的 定 标 结果 中 质量 好 的 占 比 最 低 ， 仅 为 83% ; 
流程 2 的 定 标 结果 中 质量 好 的 占 比 提高 至 96% ， 但 仍 没 有 达到 100%。 相 比 之 下 ， 由 SCAMP 参与 定 
标的 流程 3 和 流程 4 所 得 定 标 结果 中 质量 好 的 占 比 均 达 到 了 100%。 这 表明 搭配 使 用 Astrometry.net 和 
SCAMP 进行 定 标 比 仅 使 用 Astrometry.net 定 标 的 效果 要 好 。 图 11 展示 了 4 BREE to IRAR AGT [a 
均 方 根 偏差 的 分 布 。 由 图 11 可 以 看 出 ,流程 2，3 和 4 所 得 定 标 结果 的 均 方 根 偏差 分 布 相 差 不 大 ， 集 
中 程度 和 平均 偏差 均 优 于 流程 1 的 结果 。 值 得 注意 的 是 ， 流 程 3 和 4 定 标 结果 的 相差 很 小 ， 说 明 在 使 
用 SCAMP 做 扭曲 修正 前 ， 是 否 对 Astrometry.net 给 出 的 线性 定 标 结 采 进行 修正 对 最 终 定 标 结果 的 准确 
度 影 响 较 小 。 图 12 展示 了 4 种 定 标 流程 的 总 耗 时 分 布 。 图 12 KH, 流程 3 的 总 耗 时 最 少 ， 大 部 分 图 
像 在 1 s 可 以 完成 定 标 ， 这 说 明 SCAMP 计算 扭曲 修正 耗 时 极 短 ; 流程 2 的 总 耗 时 最 长 ， 一 般 需 要 S~6s 
内 才能 完成 图 像 的 定 标 ， 这 是 因为 在 流程 2 P, Astrometry.net 在 第 2 次 运行 的 过 程 中 使 用 多 个 不 同 尺 
度 的 index 文件 进行 验证 并 计算 扭曲 修正 ， 因 此 耗 时 较 长 。 
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图 7 SCAMP 修正 前 后 线性 定 标的 定性 质量 分 布 
Fig. 7 Quality distributions of the calibration results before and after SCAMP calibration 
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图 8 SCAMP 修正 前 后 线性 定 标的 均 方 根 偏差 分 布 。(a) 赤 经 方向 ; (b) HAA A 
Fig.8 Distributions of RMS deviation before and after SCAMP calibration. (a) R.A.; (b) Dec 
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图 9 solve-field 进行 线性 定 标 和 SCAMP 进行 线性 定 标 修正 的 耗 时 分 布 


Fig.9 Distributions of time cost of solve-field doing linear calibration and SCAMP doing linear calibration respectively 
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Fig. 10 Quality distributions of the calibration results from 4 different processes respectively 
讨 论 


结合 3.1 节 和 3.2 节 的 对 比 结果 以 及 两 款 软 件 定 标 原理 可 知 ，Astrometry.net 相 比 于 SCAMP 的 主 


要 优点 为 在 没有 预 完 进 行 世界 坐标 系 佑 计 的 情况 下 以 较 快速 度 对 图 像 进 行 线性 定 标 ,缺点 在 于 厂 要 进 
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行 扭 曲 修正 ， 一 次 性 定 标 虽 然 耗 时 短 ， 但 质量 较 差 ， 分 两 次 定 标 则 虽然 质量 较 好 ， 但 耗 时 较 长 ; 与 
Astrometry.net 相 比 ，SCAMP 的 主要 优点 在 于 能 够 快速 基于 给 定 的 线性 定 标 结果 进行 扭曲 修正 ， 而 对 
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给 定 线性 定 标 进行 修正 虽然 能 够 显著 提升 定 标 质量 ,但 耗 时 较 长 。 搭 配 使 用 Astrometry.net 和 SCAMP 
可 以 实现 快速 且 准 确 的 天 测 位 置 定 标 。 
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图 11 4 种 定 标 流程 所 得 结果 的 均 方 根 偏 差分 布 。(a) 赤 经 方向 ; (b) KAA A 
Fig. 11 Distributions of RMS deviation resulted from 4 different processes respectively. (a) R.A.; (b) Dec 


同时 需要 指出 的 是 ，(1) 如 果实 际 拍摄 图 像 的 扭曲 程度 比 本 文中 测试 使 用 的 ZTE 图 像 的 扭曲 程度 
严重 ， 那 么 有 可 能 导致 Astrometry.net 线性 定 标 质量 太 差 ， 从 而 使 得 SCAMP 无 法 通过 直接 扭曲 修正 获 
得 足够 准确 的 定 标 ， 此 时 很 可 能 需要 先进 行 一 次 线性 定 标的 修正 再 进行 扭曲 修正 ; (2) 本 次 实验 所 选 
择 的 测试 图 像 均 为 质量 较 好 的 图 像 ， 对 于 深度 较 浅 、 视 宁 度 较 差 的 图 像 (图 像 的 科学 价值 会 随 之 降 
IR), 不同 定 标 方式 之 间 的 对 比 结果 可 能 发 生变 化 ; (3) 本文 在 使 用 Astrometry.net 的 solve-field 命令 时 
输入 的 搜索 中 心 坐标 为 根据 图 像 已 有 世界 坐标 系 信息 计算 的 准确 中 心 坐标 ， 而 实际 数据 处 理 中 往往 只 
能 知道 图 像 大 致 的 指 癌 坐标， 因此 实际 的 解 算 耗 时 可 能 更 长 。 
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图 12 4 种 定 标 流程 总 耗 时 分 布 


Fig. 12 Distributions of time cost of 4 different processes respectively 


5 总结 与 展望 


本 文 介绍 了 Astrometry.net 和 SCAMP 两 款 软 件 的 天 测 位 置 定 标 原 理 ， 以 ZTF 的 100 张 巡 天 科学 图 
像 为 测试 图 像 ， 从 计算 耗 时 、 图 像 星 表 与 参考 星 表 匹 配 情况 、 坐 标 均 方 根 偏差 3 个 方面 比较 了 使 用 两 
牧 软 件 以 4 种 不 同 的 定 标 流程 进行 定 标的 结果 。 研 究 发 现 ， 虽 然 Astrometry.net 能 人 够 在 没有 预先 世界 
坐标 系 信 息 估 计 的 情况 下 快速 地 一 次 性 进行 线性 定 标 和 扭曲 修正 ， 但 定 标准 确 性 逊 于 耗 时 较 长 的 分 两 
步 计算 ; 而 采用 Astrometry.net 计算 线性 定 标 搭配 SCAMP 进行 扭曲 修正 的 方法 ， 可 以 做 到 在 耗 时 和 准 
确 性 上 均 优 于 仅 使 用 Astrometry.net 进行 定 标 。 在 后 续 的 研究 工作 中 ， 我 们 可 以 进一步 探究 是 否 存在 
更 好 的 选 源 方 式 、 更 好 的 软件 参数 设置 以 进一步 提高 定 标 的 速度 和 准确 性 。 
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Abstract; In order to explore the differences in the applications of Astrometry.net and SCAMP in astrometric 
calibration as well as the differences in the calibration results obtained by running them in different processes , 
we first introduce the working principles of the two suits of software, and then do astrometric calibration to 100 
images from the ZTF survey in 4 difference processes, taking Gaia DR2 catalog as the reference catalog. The 
calibration results including those given by ZTF itself are compared in terms of total time consuming, catalog 
matching, and coordinate RMS deviation. The results show that Astrometry.net can quickly calculate rough 
linear calibration results, while SCAMP can quickly perform distortion correction based on the former. The 
combined calibration process averagely takes only 1 second, and the average right ascension and declination 
RMS deviations are less than 70 milliareseconds, which is better than the calibration result only using 
Astrometry.net. Therefore, using Astrometry. net and SCAMP in combination can achieve faster and more 
accurate astrometric calibration, which may be applied to data reduction of time-domain projects nowadays. 
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